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Abstract 

For a simple one dimensional lattice version of a travelling wave equation, we 
obtain an exact relation between the initial condition and the position of the front 
at any later time. This exact relation takes the form of an inverse problem: given 
the times at which the travelling wave reaches the positions n, one can deduce 
the initial profile. We show, by means of complex analysis, that a number of known 
properties of travelling wave equations in the Fisher-KPP class can be recovered, in 
particular Bramson’s shifts of the positions. We also recover and generalize Ebert- 
van Saarloos’ corrections depending on the initial condition. 


1 Introduction 


The study of the solutions of partial differential equations describing a moving interface 
from a stable to an unstable medium is a classical subject PQ El El m [S] in mathematics, 
theoretical physics and biology P El El E]. The prototype of such equations is the 
Fisher-KPP equation (after Fisher jlUj and Kolmogorov-Petrovskii-Piskunov HU) 


du 

~di 


d‘^u 

dx"^ 


+ f{u), 


( 1 ) 


where the field u satisfies 0 < u{x, t) < 1 and where f{u) > 0. The unstable medium 
corresponds to u = 0 (i.e. /(O) = 0 and /gO) > 0) and the stable one to tt = 1 (i.e. 
/(I) = 0 and /'(I) < 0). 

One can show that equations of type 0 exhibit a continuous family of travelling 
wave solutions 

u{x,t) = W^{x — vt) (2) 

indexed by their velocities v. Explicit expressions of the travelling waves are in general 
not known except for particular velocities m- The best known example, due to Ablowitz 
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and Zeppetella [13], is w = [1 + cexp[(x — vt)/y/Q\\ for the Fisher-KPP equation Q 
with /(u) = u — and v = 5/\/6. 

Apart from describing the shapes of these travelling wave solutions a central 
question is to understand how the long time behavior of the solutions of (1^ depends on 
the initial condition u{x, 0). In general this asymptotic regime is controlled by the rate 
of the exponential decay of this initial condition. A brief review of the properties of the 
travelling wave solutions of 0 and on the way the position and the asymptotic velocity 
of the solution depend on the initial condition is given in Section 

In the present paper we study a simple one dimensional lattice version of a travelling 
wave equation. In this lattice version we associate to each lattice site n £ Z a positive 
number hn{t) which plays the role of the field u{x,t) and these hn{t) evolve according 
to 


dhnjt) _ Iah„_i(t) + h„(t) 
dt 10 


if 0 < h„(t) < 1, 

if hn(t) > 1. 


(3) 


We see that the evolution of hn{t) is linear except for the saturation at hn{t) = 1 which is 
the only non-linearity in the problem. This saturation simply means that whenever 
reaches the value 1, it keeps this value forever. The evolution Q therefore combines 
linear growth, spreading (or diffusion) because of the coupling between neighboring sites, 
and saturation, very much like in Fisher-KPP equation 0. 

The aim of this paper is to show that the evolution ^ leads to behaviors very 
similar to those expected for the usual Fisher-KPP equation 0. Moreover a number 
of properties of the solutions of Q are easier to determine than for the original Fisher- 
KPP equation (Q. Our approach is essentially based on the exact relation ( |32[ ) derived 
in Section 1^ which relates the times tn at which hn{t) reaches 1 for the first time to the 
initial condition hn{0)- We show in Section that from (32) one can obtain a precise 
description of the shape of the travelling wave solutions, in particular explicit formulas 
for their asymptotic decay. We also show in Section that ^ shares with the Fisher- 
KPP equation most of the properties expected for the dependence of the position of the 
front on the initial condition. Our results are summarized in Section E 


2 Some known properties of the Fisher-KPP class 

In this section we briefly recall some properties of the Fisher-KPP equation. 

2.1 The travelling waves 

For the Fisher-KPP equation 0 the shape Wy{x) of the travelling wave 0 satisfies an 
ordinary differential equation 

W; + uW' + /(IW) = 0 (4) 

with the boundary conditions IW(—oo) = 1 and ly^(-l-oo) = 0. By linearizing S for 
small Wy (when x is large), 

W” + uW' + f{0)Wy = 0, (5) 
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one can see that, generically, VF„(x) vanishes exponentially as x —)■ oo 

W,{x) ~ e-^", (6) 

with 7 related to the speed v of the travelling wave by 

^( 7 ) =7+—. (7) 

7 

This relation shows that depending on v, the rate 7 of the exponential decay is either 
real or complex, and these two regimes are separated by a critical velocity Vc where ^( 7 ) 
is minimum 

Vc = v{'yc) where v'i'jc) = 0 . ( 8 ) 

With v{'y) given by 0 , one gets Vc = 27 c and 7 c = a//'( 0). Under certain conditions 
on the function f{u) (such as 0 < f{u) < uf'{0) for all u see [H El [T^ and references 
therein), it is known that: 

• For 0 < V < Vc, the solutions 7 of the equation v{'y) = v are complex. The 
corresponding travelling waves solutions of (|^ oscillate around 0 while decaying 
as X —)■ 00 . 

• For V > Vc, the travelling wave is monotonically decreasing and decays for large x 
as 

VFy(x) ~ with ^ > 0, (9) 

where 71 is the smallest solution of v{'y) = v. 

• For V = Vc, the equation v{'y) = Vc has a double root 7 c and the travelling wave is 
monotonically decreasing and decays for large x as 

Wi;^(x) ~ with A > 0 . ( 10 ) 


Remark: The facts Q and (10) for v > Vc are not obvious and cannot be understood from 
the linearized equation ([^ only. These are properties of the full non-linear equation (|^, 
which can be proved under known conditions on the non-linearity f{u) (such as 0 < 
f{u) < uf'{0)). Fronts which satisfy these properties are called pulled fronts. 

For well tuned non-linearities (which fail to satisfy these conditions), travelling waves 
for V > Vc might not be monotone and the asymptotics ([^ and (10 1 might be modified; 
for instance in (§, depending on the value of v, one could have ^4 < 0 or a decay in 
exp(— 72 x) where 72 is the largest solution of ^( 7 ) = v. Rather than (10), one could 
have ^exp(— 7 cx) without the x prefactor. In all these cases, the front equation is then 
said to be pushed [El HU- 
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2.2 The selection of the velocity 


The travelling waves IT,, solutions of move at a constant speed with a time in¬ 
dependent shape. For general initial conditions u{x,0), the shape of the solution is 
time-dependent and the question of the selection of the speed is to predict the asymp¬ 
totic shape and velocity of the solution u{x,t) in the long time limit. For initial profiles 
decreasing from tt(—oo,0) = 1 to ri(-|-oo,0) = 0 it is known since the works of Bramson 
m El El ESI I under which conditions the shape of the solution u{x, t) converges to a 
travelling wave Wy solution of Q in the sense that one can find a displacement Xt such 


that 


u(Xt + x,t) Wy{x), 


with —- 
t 


V. 


( 11 ) 


In particular it is known that if the initial condition u{x, 0) satisfies for large x: 


• u{x, 0) ~ e with 0 < 7 < 7c, 

then the asymptotic velocity is ^(7), the asymptotic shape is Wyt^^-^ and 

Xt = v{^)t + Csi. (12) 


• m(x,0) <C x“e for some a < —2 (in particular for step initial conditions), 
the asymptotic velocity is Vy = i’(7c), the asymptotic shape is Wy^ and 


Xt = Vet -in t -|- Cst. 

27c 


(13) 


• ?x(x,0) ~ with a > —2, 

the asymptotic velocity Vy and shape Wy^ are the same as in the previous case but 
the logarithmic correction to the position Xt is modified: 


1 — OL 

Xt = Vyt - in t -|- Cst 

3 1 

Xt = Vyt - in t -|-in in t -|- Cst 

27c 7c 


for a > —2, 

(14) 

for a = —2. 

(15) 


(Initial conditions decaying too slowly would not lead to a travelling wave.) 

Notice that the solutions Wy of (|^ can always be translated along the x axis, so the 
“Cst” in (12 15) depends on the particular solution of that was chosen. It is often 
convenient to single out one particular solution Wy of (0: for example one may select 
the solution such that ITy(O) = 1/2 or such that / xIlT(x)dx = 0. Once a particular 
prescription for Wy is chosen, the “Cst” in the equations above is well defined. It can be 
computed in some cases such as (12), but its analytic expression is not known in some 
other cases such as (13). 
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2.3 Vanishing corrections 

The convergence property 0 does not allow to define the displacement Xt to better 
than a constant: if Xt satisfies 0, then Xt + o(l) also satisfies 0. It is however 
quite natural to choose a particular Xt, which one might call the position of the front. 
A possible choice could be 

u{Xt,t) = c, (16) 

where c G (0,1) is a fixed given number. Another possible choice would be to interpret 
—du/dx as a probability density and pick Xt as its expectation: 


Xt = — / dxx 


du 

dx 


(17) 


Either definition (16) or (17) gives a position Xt which satisfies (11). With such a 


precise definition for Xt as (16) or 0, it makes sense to try to improve on ( |12j|15| ) and 
determine higher order corrections. Ebert and van Saar loos [mils] have claimed that 
for steep enough initial conditions, the first correction to 0 is of order t and is 
universal: it depends neither on the initial condition, nor on the choice of (16) or 0, 
nor on the value c in (16), nor on the non-linearities. They found that 


Xt = Vet - 




In t -|- Cst — 3 



- 1/2 


+ 


(18) 


2.4 The Fisher-KPP class 


The main ingredients of the Fisher-KPP equation (Q which lead to travelling waves 
and fronts converging to those travelling waves are a diffusive term, a growth term and 
a saturation term. There exist many equations with the same ingredients which share 
the above properties ([8 18) of the Fisher-KPP equation: the equation satisfied by the 
travelling waves the dispersion relation 0 and the values of 7 c and Vc are modified, 
but everything else remains the same. 

To give an example which appears in the problem of directed polymers on a tree [B] , 
let us consider an evolution equation of the type 


G{x,t + 1) = j G{x + e)^p{e) de. 


(19) 


(In the directed polymers context, B is the branching ratio on the tree and p{e) is the 
distribution of the random energies associated to edges of the tree). Then u{x,t) = 
1 — G{x, t) satisfies a discrete time evolution equation with an unstable uniform solution 
u = 0 and a stable one ir = 1 as in 0 - Even though 0 is continuous in time while ( |19[ ) 
is discrete, they have similar properties: travelling waves for (19) are solutions of 


Wvix “ = y WT(a: -I- e)^p{e) de 


( 20 ) 


instead of 0. By linearizing the evolution of G{x, i) around the unstable uniform 
solution G = 1 and by looking for travelling wave solutions of this linearized equation 
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of the form 1 — G{x, t) 
replaces Q: 


exp[— 7 (x — ^( 7 )^)], one gets a new dispersion equation which 


vH) = — in 
7 


B 


=7': 


/ 9 (e) de 


( 21 ) 


but all the above behaviors (|8p(8) remain valid with Vc and 7 c computed from ( 21 ) 
and Q. For example, for B = 2 and a uniform p{e) on the unit interval (i.e. p{e) = 1 
for 0 < e < 1 and p{e) = 0 elsewhere) one gets, v{'^) = 7111 ~ 1)] which leads to 

Vc ^ 0.815172 and 7 c ~ 5.26208. 

Example (19) is a front equation where time is discrete. One could also consider 


travelling wave equations where space is discrete, say x G Z. For instance, one could 
discretize the Laplacian in Q or take (19) with a distribution p{e) concentrated on 
integer values of e. When space is discrete, special care should be taken: it is clear 
from ([^ that while the front u{x,t) lives on the lattice, the travelling wave Wv{x) is 
defined for all real values x, and even when ([^ holds, the shape of the front Wy{x — vt) 
measured on the lattice evolves periodically in time with a period \/v. Furthermore, 
the convergence (0 no longer makes any sense. One can still try to define a specific 
position of Xt by something like the following generalization of 0 : 


Xt = '^ x[u(x, t) - u{x + 1 , t)], 

x<^'L 


( 22 ) 


but with such a definition, even if the front is given by the travelling wave Wv{x — vt), 
the difference Xt — vt is no longer constant but becomes a periodic function in time 
because the shape of the front on the lattice evolves also periodically. Similarly, in the 


discrete space case, the Cst term in all the asymptotics (12 15) is in general replaced by 
a periodic function of time. 

An alternative way to locate the front when time is continuous and space is discrete 


is to invert the roles of x and t: instead of defining Xt by u{Xt, t) = c as in (161, one can 
define tx as the first time when the front at a given position x reaches a certain level c: 


u{x,tx) = c. 


(23) 


Note that when time and space are continuous, the functions Xt and tx are reciprocal 
and one can write (12][l5) as 


tx - 


X 


u(7) 


+ Cst' 


Vc 27cUc 

X ^ 1 — a 
Vc 2'ycVc 


Vc 2-/cVc 




for 

u{x, 

0 ) 

~ e-T'^ 

with 0 < 7 < 7 c, 

(24) 

In X 

+ Cst', 

for 

u{x, 

0 ) 

<C x"e" 

for some a < 

- 2 , 








(25) 

In X 

+ Cst', 

for 

u{x, 

0 ) 

~ x“e- 

for a > — 2 , 

(26) 

In X 

-^ Inlnx + Cst', 

for 

u{x, 

0 ) 

~ x“^e 

-'-fcX 

(27) 


IcVc 
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and, for steep enough initial conditions, one can write (18) as 


X 1 

tx = - 1 - 

Vc IcVc 


3 1 ^TT-Cc 

-Inx + Cst +3W 

2 w 


n-V 2 


+ 


(28) 


The main advantage of ( 24][^ | over (12 15|18 ) is that they still make sense when space 
is discrete (with a real constant Cst , not a periodic function of time). We will see that 
they remain valid for our lattice model 


3 The key formula for the position of the front 


In this section we consider the front hn{t) defined by ^ and we establish relation (32) 
between the initial condition /i„(0) and the first times tn at which hn{t) reaches the value 
1. Here we limit our discussion to the case a > 0 and to initial conditions of the form 


hn( 0 ) = 


1 for n < 0, 


^kn for n > 1, 

where the kn are non-negative, smaller than 1 and non-increasing, i.e. 


(29) 


1 > > /C2 > /C3 > • • • > 0. 


(30) 


Clearly, as o > 0, for a monotonic initial condition (30), the solution hn{t) of ^ remains 
monotonic at any later time. One can define tn as the time when hn{t) reaches 1 for the 
first time (i.e. hn{t) = 1 for t >tn while hn{t) < 1 for t < tn)- The monotonicity (30) of 
the initial condition implies the monotonicity of the times tn 


0 < ti < t2 < 


< tn < 


(31) 


Most of the properties of the solutions of ^ with the initial conditions (291 discussed 
in this paper will be based on the following exact formula 


OO A I 1 oo 


(32) 


which relates the generating function of the initial condition {kn} to the times {tn}- 


Formula (32) can be derived as follows. If one defines the generating functions 

Hm{t) = hn{t)X^-^+\ 


(33) 


n>m 


one can see that for tm-i < t < tm (with the convention that to = 0) the evolution of 
Hmit) is given by 

= (1 + a\)Hm{t) + a\. (34) 
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This of course can be easily solved to give 

Emit) = -rr^ + (35) 

1 + oA 

where the $m’s are constants of integration. These can be determined by matching 
the solutions at times 0, ti, t 2 , ■ ■ ■ '■ 

OO 

El{0) = ^ knX"', = A(l + ■, (36) 

n=l 

and one gets that for tm-i < t < tm 


Emit) = - 


aX 


I 1 OO 

n=m 


1 -|- clX 1 -|- dX 


(37) 


Then (32) follows by taking m = 1 and t = 0 in 


Remark that formula (32) appears as the solution of a kind of inverse problem: given 


the times tn, one can compute the initial profile kn by expanding in powers of A. This 
gives expressions of kn in terms of the times tm's for m < n. Alternatively one can 
determine the times tn in terms of the initial profile km for m < n: 

^ I k\ /'I ~ \ ^ 1^^ \ ^ \ 


e -^1 = 


a + 1 ’ 


aki + k2 , ^ 

e ^ ,-h atie 


g ts = 




Cl -|- 1 


(38) 

etc. Unfortunately these expressions become quickly too complicated to allow to deter¬ 
mine how the times tn depend asymptotically on the initial profile {kn} for large n. How 
these asymptotics can be understood from (32) will be discussed in Section 


4 Travelling wave solutions 


4.1 The exact shape of the travelling waves 


As usual, with travelling wave equations, the first solutions one can try to determine are 
travelling wave solutions moving at a certain velocity v. Because the hn{t) are defined 
on a lattice, a travelling wave solution moving at velocity v satisfies 

hn{t) = hn+i(t + (39) 


Clearly this implies that the times tn form an arithmetic progression, and by shifting 
the origin of time, one can choose 


n 

tn — • 


(40) 


This immediately gives, using (371, the generating function of the front shape at all 
times: for example for 0<t<ti = l/u, one takes m = 1 in (37) and gets 


E ^n{t)X^ = 

n>l 


aX a + 1 

1 -|- aA 1 -|- aA ^ g(i+aA)/'u _ ^' 


(41) 
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Another way of determining the travelling wave solutions is to look directly for 
solutions of Q of the form (391. One sees that Wv must satisfy 


Wv{x) = 1 for X < 0, W^{x) + aWy{x — 1) + vWy(x) = 0 for x > 0. 


(42) 


These equations can be solved iteratively: for x < 0, one already knows that Wy{x) = 1. 
For X G [0,1] one has therefore Wy + vWy + a = 0, which implies for (x G [0,1]) that 
Wy{x) = (a + — a (the integration constant being fixed by continuity at x = 0). 


Knowing Wy{x) for x G [0,1], one can solve (42) for x G [1,2] and so on. 


Wy{x) = 


(a + l)e 


—xjv _ 


if X < 0, 
if X G [0,1], 


a(i+a ) ^ _ ^^g(i a;)A _|_ (X + a)e ifxG[l,2], 


(43) 


In fact one can solve directly (42| by considering for x G [0,1] the generating function 


(A, x) I—A"lFn(x + n). Then, as can be checked directly from (411 and (42), Wv{x) 
and hn{t) are related for all n G Z and t > 0 by 


hn{t) = Wy{n - vt). 

4.2 The decay of the travelling waves 


(44) 


The large n behavior of the travelling wave hn{t) (or equivalently the behavior of Wy{x) 
for large x) can be understood by analyzing the singularities in A of the right hand side 


of (411. These singularities are poles located at all the real or complex zeros of 

g(l+aA)/,; - A = 0. 


(45) 


(one checks there is no pole at A = —1/a) and each pole gives rise to an exponential 
decay in hn{t). Using A = exp( 7 ), ([45|) can be rewritten into 


vi'j) = 


1 + ae"^ 
7 


(46) 


which is the dispersion relation for (|^, similar to ([^ or to (21). In fact, one can obtain 
(46) as in Section]^ by looking for the velocity v compatible in (42) with an exponentially 
decaying travelling wave of the form Wy{x) ~ 

One is then led to distinguish three cases depending on the number of real solutions 


of (46). There is a critical value Vc where (46) has a double zero on the real axis. This 


critical velocity Vc and the corresponding decay rate jc are the solution of 


a(7c - l)e^= = 1, 


Vr. = 


1 

7c - 1' 


(47) 
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1. For V < Vc, there is no real A solution of (45), but there are complex roots. The 


large n behavior of hn{t) is then governed by the two roots Ai and A^ of (451 closest 
to the origin 

Because Ai and A5( are complex, hn{t) changes its sign as n varies. So for v < Vc, 
as in the Fisher-KPP case, there are travelling wave solutions, but they fail to give 
positive profiles. 


2. For V = Vc given by (47), there is a double real root Ac = of (45). Then for 

(49) 


large n the profile is of the form 

2(1 +a) 


hn(0) ~ 


I+ Vc 


n 


1 + 4i;c 


3(l + r;c)J 




3. For V > Vc there are two real roots 1 < Ai < A 2 of (45) and the large n behavior 
is controlled by the smallest root: 


hniO) ^ 


(a + l)v 


(1 + aAi)(t' — aAi) 


Ar” 


(50) 


We see that for all velocities, we get explicit expressions of the prefactors of the expo¬ 
nential decay of the travelling waves. These prefactors are in general not known for more 
traditional travelling wave equations, such as the Fisher-KPP equation ([^. 

In each case, corrections to (48 50) can be obtained from the contributions of the 


other roots of (45). For instance, in the cases v < Vc oi v > Vc one could write 


^n(O) cr; 


(a + l)t; 


(1 + a\r){v — aXr) 


Xr 


(51) 


where the sum is over all the complex roots Ar of (451. In Figure we compare the exact 
solution (43) of (42) with the asymptotic expansion truncated to a finite number of 
roots of (451 closest to the origin and one can see that the truncation gives a very good 


fit of the actual solution. 

We have seen that the travelling waves for v < Vc were oscillatory. For v > Vc, they 


decrease monotonically towards 0; this can be seen directly from equation (42) verified 


by W^{x)\ write Wy{x) = R{x)e with 7 a real positive number related to v through 


Then (42) gives 


the dispersion relation (46). (Notice that choosing such a 7 is impossible if f < Vc-) 


R{x) = for x < 0, ae'^[R{x — 1) — R{x)] + vR'{x) = 0 for x > 0. 


(52) 


As R{x) is strictly increasing for x < 0 it must be strictly increasing for all reals; 
otherwise, on the first local maximum Xm, one would have R'{xm) = 0 and R{xm) > 


R{xm, — 1) which is incompatible with (52). Hence, Wy(x} is positive and, from (42) 
strictly decreasing. 


10 



























Figure 1: The travelling wave Wv{x) solution of (42| for v = 4 and a = 1 as a function 
of X. The plain line labeled “exact” is the exact small-x solution (43). The dashed lines 
are the sums (511 truncated to a given number of first terms: “with 2 roots” means only 
the two real roots Ai and A 2 , “with 4 roots” means the two real roots and the first pair 
of complex conjugate roots and “with 22 roots” means the two real roots and the ten 
pairs of complex conjugates roots closest to the origin. The inset is a zoom of the small 
rectangle around x = 0 and Wy = 1. 


5 How the initial condition determines the asymptotic regime 


We now discuss how the position of the front at large times (or equivalently the large n 
asymptotics of the times tn) depends on the initial condition. 

First, by using mostly a comparison property, we will show that the final velocity of 
the front is determined by the large n decay of the initial condition kn- Then, we will 
recover the logarithmic corrections ( I2|[l5 ) and sub-leading terms as in (18) by analyzing 
the key relation (32) between the initial profile kn and the times tn- 


We write (321 as 


(1 + aX)K{X) = -aX + {a + 1)T{X), 
where the two functions K{X) and T{X) are defined by 


(53) 


K{X) = £ knX^, r(A) = £ e-(l+«A)in^n_ 


(54) 


n=l 


71=1 


The large n behavior of the knS and of the determines the domain of convergence 


of these two sums and one can try to use (53) to relate their singularities. 


When X = > 1, we will often use the following form of T(X) written in terms of 

the dispersion relation v{'y): 


r(e^) = (55) 

n=l 
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5.1 Selection of the velocity 

Let us first show that the final velocity of the front is determined by the large n behavior 
of the initial condition in the same way as for other equations in the Fisher-KPP class. 
To do this, we use an obvious comparison property; considering two initial conditions 
{kn'^} and {kn^} with the corresponding times and one has 

if 0 < k^'^ < k^'^ for all n, then for all n. (56) 

To keep the discussion simple, we focus only on initial conditions {kn} with > 0 and 
the following simple asymptotics: 

• If /cn ~ with 0 < 7 < 7 c. 

Pick an e > 0 small enough so that 0 < 7 — e and 7 + e < 7 c, and consider the two 
travelling waves going at velocities v{'y — e) and u (7 + e) (they decay respectively 
like and 6 “^'^"’“^^"'). It is clear that the initial condition {kn} can be 

sandwiched between these two travelling waves suitably shifted in space, so that, 
by using the comparison property one gets 

^ ^ < liminf — < limsup — < ^ (57) 

v{j — e) n^oo fi n^oo n v{j + e) 

Now take e —)■ 0 to get 

lim — = . (58) 

rn-oo n u( 7 ) 

• kn = 0 . 

It takes a time tn to have hn{t) = 1. But at time tn, the hn+mif) for m > 0 are 
positive so that, from the comparison property, one has 


tn+m ^ tn + tr, 


(59) 


The sequence {tn} is sub-additive and therefore tn/n has a limit which we call 1/v. 
By comparing the initial profile kn = 0 to the travelling wave going at velocity Vc, 
one must have 1/v > 1/vc- 

We are now going to show that 1/v cannot be strictly greater than 1/vc- Indeed, 
if we had 1/v > 1/vc, the series ( |55[ ) defining T(A) would be uniformly convergent 
on the whole positive real axis A because v{'y)tn/n would eventually be larger than 
1 -|- e for some e > 0. One would then get 


r'(A) = Y, 


n 


atr. 


for all A > 0 


n>l 


(60) 


because the series (60) would also be uniformly convergent. 


However, for real and large enough A (at least for A > min„ n/(at„)), one would 
obtain T'{X) < 0. But, with = 0 one has K{X) = 0 and from T'(A) = 
a/{a -|- 1), in contradiction with T'{X) < 0. 
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We conclude that one must have 

lim ^ = Jl if = 0 for n > 1 . (61) 

n-s>oo n Vc 

• If fen ~ or if kn = 

Again, by the comparison property, the initial condition can be sandwiched be¬ 
tween fen = 0 and, for any e > 0 , the suitably shifted travelling wave going at 
velocity v{'yc — e). This leads to conclude that 

lim ^ ( 52 ) 

n^oo fi 

The velocity selection thus works as for other equations of the Fisher-KPP class. 


5.2 Sub-leading corrections 

We limit the discussion to initial conditions similar to those discussed in the previous 
section which lead to a front with some asymptotic velocity V: 


,. tn 
n—j-cxD 77, 


1 

V' 


(63) 


We also assume that fen > 0 which implies that V > Vc as was shown in the previous 
section. 

If P > Vci write V = ^(71) = ^(72) with 0 < 71 < 7c < 72. We have seen that 
this velocity is reached for initial conditions such as fen ~ In (55), it is then 

clear that the series T(\) is divergent for A E (e'^Te'^^) and convergent for A < or 
A > Furthermore, in (54|, the radius of convergence of K{X) is and, as fen > 0, 
the function K(X) must have a singularity at A = We thus see that both T{X) and 
K{X) become singular as A approaches from below on the real axis. By matching 
the singularities of these two functions, we will obtain the sub-leading corrections to 
for large n. 

For V = Vc, if the initial condition is kn ~ the same argument applies: both 

r(A) and K(X) are singular when A reaches and one must match the singularities. 
But, with V = Vc, one could also have an initial condition which decays faster than 
and for which the radius of convergence is larger than (even, possibly, infinite). Then, 
of course, K(X) would have no singularity at A = even though the convergence of 
T(A) would remain problematic when A approaches We will see that the large n 
behavior of tn is tuned to “erase” the singularities in T{X) at to satisfy (53). 

We attack the problem by assuming that the tn are given and we try to obtain the 
asymptotics of the kn- The starting point is thus to assume a velocity V = ^( 71 ) with 
7 i < 7 c) and study T{X) when A gets close to e'>'T If one chooses, in all generality. 


t 


n 


_L 

V ^ 7 iP’ 


(64) 
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where bnjn —)■ 0, one gets from (55) 


T(e") = E 




(65) 


n=l 


Now we want to take 7 = 71 — e and expand for small e in order to extract the nature 
of the singularity. Two cases arise: 

• liV > Vc (which means 71 < 7c), then f'(7i) < 0 and to leading order 

r(e'^i“^) = E exp --for y > Uc, (66) 

n=\ 

with /X = 1/71 + v'{'^i)/V. 

• l{V = Vc (which means 71 = 7c), then u^(7c) = 0 and one must push the expansion 
further: 


CXJ 

= E®^P ( 


n=l 


r/ lcv"{lc)^2. 

e + -- 


Lv 2vc 

^ ^ Ic 

/ J 


for V = Vc- (67) 


It is already clear that cases V > Vc and V = Vc need to be discussed separately. 


Equations ( 66 ) and (67) are the starting point of our analysis which is presented in detail 


in the following subsections. 

We will make heavy use of the following formulas: for e > 0 small. 


E 


n>l 


singular 


r(i + a) 

gl+O 

(-l)"e-"-Mne 


if a is not a negative integer. 


( 68 ) 


(-a - 1 )! 

r-r(l + a) lne + 0(l) 


E(lnn)n"e 


n>l 


ci+CK 


if a is a negative integer. 

if a is not a negative integer, 


= 


singular 


(- 6 ) 


—a —1 


(-a - 1 )! 


InW 


+ 0 (lne) 


if a is a negative integer, 


(69) 


where the meaning of “singular” for a function F(e) with a singularity at 0 is that the 
difference between F(e) and F(e) is a regular function of e which can be expanded 

as a power series. 


5.2.1 For V > Vc 

As explained above we write V = u( 7 i) with 71 < 7 c, and we choose such that tn/n 
that converges to 1/F. If one chooses 


n Blnn + C 


(70) 
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by keeping the leading order in ( 66 ) and using (681 one gets for B not a positive integer 


singular V— 


l-B 


(71) 


It is then easy to check that matching the singularities leads to an initial condition 
decaying as 

^ (l + a)e-^ r Vn ^ 72 ) 

Remarks: As can be easily checked, even though is not valid if i? is a positive integer 


(72) is. One can also check that for B = C = 0 one recovers the asymptotics (50) of the 


travelling wave. 


5.2.2 For V = Vc 

if F = Vc, the main difference with the previous case is that u( 7 c —e) — u( 7 c) ~ as e —)• 0 
and one must use the expansion (67) instead of ( 66 ). As before, we choose a specific 


form for the times tn which allow to easily make the comparison with the different cases 


(13 15) of the Fisher-KPP equation: 


n 

tn — 

Vc 


B\nn + C 

IcVc 


(73) 


With 5n = B\n+C into (67), one obtains generically (when B 0 {1, 3/2,2,5/2,3,...}, 
see discussion below) 


r(e^^-") 


-c 


singular 


T{l-B) 


7c^"(7c) 

2vc 


B-l 


2B-2 


(74) 


Then, using ( 68 ) again and (53), one gets 


u ^ + 1 -c ^7W"(7c) 

- t- 


IcVc 


2vc 


B-l 


r{i-B) 
T{2-2B) 


2Bg-7cn_ 


(75) 


We see that the asymptotics of the initial condition (75) and of the times (73) for large n 


are related as in the Fisher-KPP case (26) and that the constant term in (26) can be 


determined. As in (70), one must be careful when B is a positive integer: (jM) should be 
modified to include the logarithmic correction of ( 68 ), but (175^ is not modified as can 


easily be checked (the ratio of the two Gamma functions has a limit). 

There is another difficulty when B G {3/2,5/2, 7/2,...}: for these values, the ratio 
of Gamma functions in ( |75[ ) is zero. This means that an initial condition {kn} leading 
to (73) with = 3/2 (for instance) must decrease faster than For these special 

values of B, the right hand side of (74) is actually regular as g 2 B -2 jg non-negative 


integer power of e; any singular part of must come from higher order terms. 

We are now going to show that no non-negative initial condition {kn} can lead to a 
time sequence {t„} with an asymptotic expansion starting as in (73) with B > 3/2. To 
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do so, we will show that the initial condition = 0 leads to (73) with B = 3j2 (plus 
higher order corrections). As any non-negative initial condition must lead to times {tn} 
which are smaller than the times of the = 0 initial condition, this will prove that B 
cannot be larger than 3/2. 

Consider therefore the case kn = 0; one has K{X) = 0 and, from (53), one gets 
T(A) = a\/{a + 1). Obviously, T(A) has no singularity as A approaches so the right 


hand side of (74) must be regular, which implies that B G {3/2,5/2, 7/2,...}. We now 


rule out any value other than B = 3/2 by looking at the term of order e in the expansion 
of One can check that the only terms of order e come from the ebnj^c in (67) 

and from \i B = 3/2, so one has 

1/2 


= r(e^^) + 


1 


oo 


- V<5„e-^"+e-^r(-l/2) 


-Ct 


n=l 


f lcV''{lc) 
V 2Vr. 


1b=3/2 


e-l-o(e). (76) 


Notice (64) that <5^ > 0 for the = 0 initial condition because it is below the travelling 
wave at velocity Vc for which 5n = 0. The first term of order e in (76) is therefore 
positive; on the other hand, the second term (only B = 3/2) is negative. But, from 
T(A) = a\/{a + 1) the term of order e must be negative; therefore one must have 
B = 3/2 for the zero initial condition and, therefore, B <3/2 for any non-negative 
initial condition. 


To summarize, the relationship between the times (73) and the initial condition (75) 


we have established in this section is valid only for B < 3/2 because we only consider 
non-negative initial conditions. Furthermore, to have (73) with B = 3/2, one must 


have an initial condition decreasing faster than the n suggested by (|75|). No 


non-negative initial condition can lead to (73) with B > 3/2. 


5.2.3 For V = Vc and B = 3/2 


The case 7? = 3/2 is of course the most delicate and it corresponds to ( |13|15|18l ) in the 
Fisher-KPP case. For the tn given by (73) the leading singularity is not (74) but rather 


r(eT'=-^) 


singular 


V IcVc 


(77) 


(This term comes from the first order expansion of the ebn/'^c term in (67).) Relating 


this to the {kn} through (53), it leads through (68) to ~ n with a negative 


prefactor. So there is no way for a non-negative initial condition to be compatible with 


exactly (73), without any extra term. 


Therefore, we need to add some corrections to (73) when B = 3/2. Let us consider 
a correction of the form 

n I In n -|- C -|- Dn~^ 


tn — -h 

Vr. 


IcVc 


for some ^ > 0. Plugging this correction into (67) one gets 

3e 




(7c) _3 

3 2uc n 2 


n=l 


1 -|- In n — Dn ^ -|- • • • 
27c 


(78) 


(79) 
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where the “•••” contains smaller order terms of orders ne^, en“^, e^ln^n, etc. 

Consider in turns the terms in the square bracket. The “1” leads to the right hand side 
of (74) with B = 3/2, which is simply a regular term linear in e. The term in elnn gives 
the right hand side of (77) and the —Dn~^ contribution can be computed from 


oo 

Ee- 

ri.=l 


-ycV (7c) ^2^ 3 c 

2vc ^ ^ n~2~^ 


singular 


il ( = I 

(80) 

Several subcases must be considered 

• If ^ > 1/2 this is smaller than In e; therefore the leading singularity is still given 
by ([^ which is incompatible with a non-negative initial condition. 

• If 0 < ^ < 1/2 the leading singularity for is as given by (80). This 

leads to 


kn — —De 


_c r(-i-g) / 7 c^"(7c ) ^ 

7cUc r(-l - 2^) V 2uc 


^-2-2^g-n7c_ 


(81) 


With 0 < .^ < 1, this is positive if H > 0. 


If ^ = 1/2 the corrections from (80) and from (77) are both of order e^lne. This 
leads to 


\ Qj _ ^ 

kn ^ 2 - e ^ 

IcVc 


dWM_3, 

Vc ' 


/27ru"(7c 

IcVc 




(82) 


which is positive if D is large enough. Notice also that the square bracket in (82) 

vanishes for ^_ 

27rUr 


^ = 3W 3 y (83) 

V Ty\ic) 

This means that initial conditions decaying faster than (including the 

zero initial condition) must lead to (78) with ^ = 1/2 and D given by (83). This 
is exactly the prediction (28). 

To finish, notice that we found the first terms of the asymptotic expansion for the 
times tn when the initial condition decays as when a > —2 (see (73) for B < 3/2) 

and when a < —2 (it is of the form (78) with ^ = —1 — a/2 for —3 < a < —2 and ^ = 1/2 
for a < —3), but we did not yet considered the case where kn — n~‘^. One can 
check that by taking, as in (27), 


n 


tn — -1" 


In n — In In n -|- C 
IcVc 


one obtains 


r(e^=-") 


singular 


~ e 


-C /8vr7cu"(7c) 


eln e. 


(84) 


(85) 
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which leads to 


87r7c?;"(7c) 


IcVc 




(86) 


6 Summary 


In the previous section, we have computed the initial conditions kn as a function of the 
times tn- Table summarizes our results. 


1 

kn ~ n°‘e 
with 7 < 7c 

71 1 r 1 

~ , , -I- , , a\n.n + C 

vin) 7v{j) L J 

see i 

o 

72 

1 

2 

kn ~ 

with a > —2 

n 1 

f -1- 

Vc IcVc 

' 1 — a , ' 

-In n -|- C 

2 


see 

73 

75 


3 

kn ~ 

n 1 

f -1- 

Vc IcVc 

3 1 

- In n — In In n -|- C 

[2 J 

see 

84 

86 

1 

4 

kn ~ 

with —3 < a < —2 

n 1 

f -1- 

Vc IcVc 

"3 a~\ 

- In n -|- C -|- 

[2 J 

see ( 

78 

81 

82 

1 

5 

kn <€. 

for some a < — 3 

n 1 

tn^ + 

Vc IcVc 

— In7i-|-C* 3^ 

2'kVc _i 

lcv"{'yc) ^ 


see ( 

78 

82 

83 

1 


Table 1: Asymptotics of tn as a function of the leading behavior of the initial condi¬ 
tion kn- 


These asymptotics agree with all previously known results discussed in Section 
Case 4 is a new prediction, and the domain of validity of Ebert-van Saarloos correc¬ 
tion (18) from m is made precise (case 5). 

The constant C can easily be computed in cases 1 to 3, but we did not manage to get 
a closed expression in cases 4 and 5. Similarly, we have no expression for D in case 4; in 
particular, for a = —3, the D of case 4 is not given by the prefactor of in case 5 

because for a = — 3 the right hand side of (82) must not vanish. 

The vanishing terms and in cases 4 and 5 depend only on the leading 

behavior of kn for large n. One could compute higher order corrections in cases 1 to 3 
using the same technique by looking at the next singularities in T(e'^“'^), but one would 
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need to know a bit more about the asymptotic behavior of kn- one would find that 


If kn 



in case 1 for a / 0, 

in case 2 for a 0 { — 1,0,1}, (87) 

in case 3, 


where the prefactor D could be computed in each case. These vanishing corrections in 
cases 1 to 3 are less universal than in case 4 to 5 as they do not depend only on the 
leading behavior of kn for large n, but also on the fact that the sub-leading behavior 
of kn decays fast enough compared to the leading behavior. For case 1 with a = 0 and 
case 2 with a = 1, the initial condition behaves asymptotically as the travelling wave 
eventually reached by the front, and vanishing corrections might depend on the initial 
condition in a more complicated way. Case 2 with q: = —lora = 0 are border cases 
with slightly different corrections. 

If we conjecture that the new results (cases 4 and 5) of Table hold for the whole 
Fisher-KPP class one can obtain, by inverting the relations between tn and n of Table 
the asymptotics of the position Xf for initial conditions of the form u{x,0) ~ 

This is done in Table [H 


1 

u{x, 0) ~ x"e 
with 7 < 7c 

Q/ 

Xt ~ u( 7 )t H— In t -|- C' 

7 

2 

u(x, 0) ~ 
with a > —2 

(y. — 1 

Xt ~ Vet + — - Int + C' 

3 

tt(x, 0) ~ x~‘^e~^'"^ 

3 1 

Xt — Vet -In t H-In In t 

27c 7c 

4 

u(x,0) ~ 

with —3 < a < —2 

Xtc^Vct- — Int + C' - 
27c 

5 

u{x,0) <C 

for some a < — 3 

3 , / 27r 1 

XtC:tVct \D.t + C 3J t 2 

27c V 7c^^ (7c) 


Table 2: Asymptotic expansion of Xt as a function of the leading behavior of the initial 
condition m(x, 0). 
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7 Conclusion 


The main result of the present work is the exact relation (321 between the initial condition 


and the positions of the front at time t for the model (3). Relating the asymptotics of 
the tn's to those of the fc„’s, using the exact relation (32) is an interesting but not 


easy problem of complex analysis. It allows to obtain precise expressions of the shape 
of the travelling waves, including prefactors which are usually not known in the usual 
equations of the Fisher-KPP type. It also allows one to recover the known long time 
asymptotics of the front position, and to get previously unknown results; in particular, we 
have shown how fast an initial condition should decay to exhibit the Ebert-van Saarloos 
correction, and that there is a range of initial conditions which exhibit the —3/2Inf 
Bramson logarithmic term but for which the Ebert-van Saarloos correction is modified. 
(See cases 4 and 5 of Tables and H) 

As shown here the analysis of the asymptotics using complex analysis, is tedious 
but rather straightforward. Higher corrections to the asymptotics of the position could 
be determined. One could also try to study how, depending on the initial condition, the 
asymptotic shape is reached. Furthermore, it would be interesting to generalize (© to 
evolutions involving more than two neighboring sites, or to a non-lattice version of the 
model. More challenging would be to attack the noisy version of the problem [IHIEDI. 


A An heuristic derivation of the positions of the front 

In this appendix we show that several expressions of the position of the front for a 
Fisher-KPP front can be recovered by considering a simplified version of the Fisher-KPP 
equation (Q where the non-linear term is replaced by an absorbing boundary. Consider 
the following linearized Fisher-KPP equation with a given time-dependent boundary Xt 
with Xq = 0: 

A = g + /(o)„ 

u{Xt,t) = 0. 

For a given y > 0, we look at the value u{Xt + y, t) of the solution at a distance y from 
the boundary. Intuitively, if Xt increases too quickly with t, this quantity is pushed to 
zero. On the other hand, if Xt increases too slowly, it diverges with t. It is only for 
finely tuned choices of Xt that u{Xt + y, t) remains of order 1. 

Now we suppose that Xt is no longer given a priori but is instead determined by 


u{Xt + l,t) = 1. 


(89) 


It has been shown [2T] that the solution of ( 88|89 ) for compactly supported initial con¬ 
ditions leads to the same long time asymptotics for Xt as for the Fisher-KPP equation 
(see Section]^: one recovers the Bramson term (13) and the Ebert-Van Saarloos correc¬ 
tion ( [T^ . 

For initial conditions decaying fast enough, one expects Xt to be asymptotically 
linear. If Xt were really linear (not only asymptotically but at all times), (88) would 


20 











be very easy to solve. In this Appendix, we solve a simplified version of ( 88 ) where the 
boundary is replaced by a straight line. This allows us to recover the velocity and the 
logarithmic corrections ( 12](l5 ) of the Fisher-KPP equation. 

The version we actually consider is therefore the following: For each given time t, we 
replace the boundary by a linear boundary of slope Xt/t and solve 


du d'^u , V 

a^ = a^ + /(o)“ 

t 


u{^s,sj = 0 . 


(90) 


We then tune the value of Xt to satisfy (891 at time t. 
m initial con 

g{x,s\xo) = 


For an initial condition 6{x — xq) the solution to (90) is 
eT(o)s 


exp 


(x - Xo)" 


4s 


exp 


Xt {x + xqY 


t 


■xo 


Taking s = t and writing x = Xt + y, one obtains 


+ y,t\xo) = 


1 


\/ATrt 


exp 


rmt - 


(Xt + yY + xl- 2XtXo 
4t 


4s 


2 sinh 


V 2t 


Given a general initial condition u(xo,0) for xq > 0 one has 

roo 

u{Xt + y,t) = / dxog{Xt + y,t\xo)u{xo,0), 
Jo 

which, after writing Xt = ct — 5t with 5t <C t, leads to 


(91) 


(92) 


(93) 


u{Xt + y,t) = 


y/Yi 


exp 


with It{y) = / dxo m(xo, 0 ) exp 


(y-^t 


(y-S 

} 

4t 

t^txo 

2t 

Xq 

4t 

sinh 1 


X h{y), 


(94) 


\ 2t r 


Depending on the initial condition ri(xo,0), we can now determine for which values 
of c and 6t the front u{Xt + y, t) remains of order 1 for y of order 1 as t increases. 

• For u(xo,0) ~ with 7 < c/2, 

one finds that the integral It{y) is dominated by xq — (c — 2^)t. One obtains 

^ r.2 


—exp + 


c — 27 


and u{Xt + y, t) ~ 2A sinh 


c — 27 


-y 


4 

exp 


(f'(0) - 7C + + 7<^t - 2^/ 


(95) 


Writing u(Xt + y,t) ~ 1 leads to c = 7 + /^( 0)/7 = ^( 7 ) and ~ Cst. The starting 
hypothesis 7 < c/2 then translates into 7 < 7 c = \//'(0). We conclude that 


For u(xo, 0) ~ e with 7 < 7 c, Xt ~ v( 7 )t + C, 


(96) 


as in ( 12 ). 
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For rt(xo,0) ~ AxQe~'^^° with 7 < c/2, 

the integral hiy) is again dominated by xq — (c — 2'y)t. The large t expression 
of u{Xt + y,t) has an extra term [(c — 27)t]" which is canceled by taking now 
5t — — - Int + Cst. (The value of c remains the same.) We conclude that 


a . 


For u(xo,0) ~ Xq e with 7 < 7c, W ~ v(j)t H— Int + C. 

7 

For rt(xo,0) <C 6“'’'^° for some 7 > c/2 (steep initial condition), 
the integral It{y) is dominated by xq of order 1. This leads to 


(97) 


y roo CXq 

h{v) - Yt Jq ’ 


and u{Xt + y, t) ~ exp 


i(/(o) - j) - x(y-<5t) 


(98) 


One needs to take c = 2yJ /'(O) = Vc = 27c and 6t = ^ Int + Cst. The starting 
hypothesis 7 > c/2 translates into 7 > 7 c and we conclude that 

3 

For tt(xo, 0) <C for some 7 > 7c, W ~ Vct -In t + C, (99) 

27c 


as in (13). 

For u(xo,0) ~ 74xoe“2*o, 

depending on the value of a, the integral It{y) is dominated by values of xq of 
order 1 or of order \/t. In any case, xq ^ t and one can simplify It{y) into 


y 

7t(y) - 7^ / dxou(xo,0)xo 
2t Jo 


exp 




( 100 ) 


When a < — 2, this integral is dominated by xq of order 1, the Gaussian term can 
be dropped and one recovers (98) and (991. 

When a > —2, the integral is dominated by xq of order \/t. One gets 

^7ly2"F(l + ifa>- 2 . 


h{y) ^ dxo Xo+^ exp 


1 

1 

loto 

J 

At 

“I 


. In t 

Ay— 
" At 


if a = —2. 


( 101 ) 


Into (94) one must therefore take c = 2 ^/f^0) = Vc = 27c and 5t = 77^ In t + Cst 
if a > —2 or 6t = nir Int — In Inf if a = —2. We conclude that 


27c 


7c 


Vet -In t + C 

27c 


if a < —2, 


Q 1T1 1T1 / 

For u(xo, 0) ~ Xoe“'^‘'^°, Xt ■{ Vet -IntH- \-C if a = —2, (102) 

27c 7c 


Vet - — In t + C 

27c 


if a > —2, 


as in (13 15). 
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